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, We describe the use of a near-field scanning microwave microscope to image the permittivity and 

\^ • tunability of bulk and thin film dielectric samples on a length scale of about 1 /im. The microscope is 

^SJ ' sensitive to the linear permittivity, as well as to nonlinear dielectric terms, which can be measured 

as a function of an applied electric field. We introduce a versatile finite element model for the 
system, which allows quantitative results to be obtained. We demonstrate use of the microscope 
O \ at 7.2 GHz with a 370 nm thick Bao.6Sro.4Ti03 thin film on a LaAlOa substrate. This technique 

. is nondestructive and has broadband (0.1-50 GHz) capability. The sensitivity of the microscope to 

changes in permittivity is Ae^ = 2 at e,. = 500, while the nonlinear dielectric tunability sensitivity 
is Aeii3 = 10"^ {kV/cmy^. 



o 



> 



I 



I. INTRODUCTION 



Dielectric thin film research has become increasingly important as the demand grows for smaller, faster, and more 
reliable electronics. For example, high permittivity thin films are under study in order to fabricate smaller capacitors 
while minimizing leakage. Low permittivity materials are being sought to allow smaller scale circuits while minimizing 
undesirable stray capacitance between wires. Nonlinear dielectrics, which have a dielectric permittivity which is a 
function of electric field, are being used in tunable devices, particularly at microwave frequencies. Finally, ferroelectric 
thin films are a solution for large-scale, non-volatile memories, 
^-j.' ' All of these dielectric thin film technologies demand high-quality, homogeneous films. However, thin film properties 
, and homogeneity can be difficult to measure. Consequently, it is important to have a set of reliable techniques for 
' evaluating thin film properties, such as permittivity arui-|tunability. A jjiumber of different techniques are available. 
, Thin film capacitors which are sensitive to the normalau and in-planecl components of the permittivity tensor can 
• be used at low frequencies 1 MHz). Microwave measurements of permittivity have been made usins^transmission 
measurements through a microstrip structurecl, and reflection measurements from a Corbino probeOo and similar 



geometries.Q However, these techniques average over large areas, and involve depositing thin film electrodes, which 

3iJj re 



itself can alter the properties of the sample. Diolectridl and opeiJj resonators |liave been used as well, but suffer 
from low spatial resolution. Cavity perturbationll2l and waveguide transmission^ techniques have also beep_.used, 
■ but have the disadvantage of measuring the entire sample. More recently, near-field microscopy techniquesE3t3 have 
Q ^ allowed quantitative measurements with spatial resolutions much less than the wavelength. These techniques use a 
O • resonator which is coupled to a localized region of the sample through a small probe, and have the advantage of being 
^ non-destructive. However, it is still difficult to arrive at quantitative results and maintain high spatial resolution. 
• i Previously, we introduced our solution, a near-field scanning microwave pajpcpscopc (NSMM), and its apftlication 
' to the low-resolution (100-500 /im) quantitative imaging of sheet resistancell3't£l and dielectric permittivity.E3 In this 
^ . paper, we present a methoji-for using the NSMM for high-resolution {^-^ 1 /im), quantitative imaging of thin film 
- - ' permittivity and tunability.lij 
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II. DESCRIPTION OF THE MICROSCOPE 



Our NSMM systems (see Fig. ^ consist of a coaxial transmission line resonator which is terminated at one end 
with an open-ended coaxial probe, and at the other end is coupled to a microwave source with a coupling capacitor. 
The probe, which has a sharp-tipped center conductor extending beyond the outer conductor, is held fixed while the 
sample is raster scanned beneath the probe. A sample holder gently presses thje sample against the probe tip with a 
force (about 50 /iNcEl) which is small enough that the sample is not scratched.E3 Because the rf fields are concentrated 
at the probe tip, the resonant frequencw-|Md quality factor Q of the resonator are a function of the sample properties 
near the probe tip. A feedback circuitll3'E3 keeps the microwave source locked onto a selected resonant frequency of 
the microscope resonator, and outputs the frequency shift (A/) and Q of the resonator. This data is stored by a 
computer, which also controls the scanning of the sample beneath the prob e ,,. W e have shown previously that the 
spatial resolutiop-oLthe microscope in this mode of operation is about 1 /im.tlHl^ 

Our probesllZHlj'LiJ are constructed using a short section of semi-rigid coaxial cable, where we have replaced the 
center conductor with a capillary tube of. the same outer diameter. The probe tip is a commercial scanning tunneling 
microscope (STM) conical tungsten tiptj which is inserted into the capillary tube, and held in place by friction. 

We apply a local dc electric field to thin film samples with a bias teeE3 in the microscope resonator (see inset to 
Fig. U), allowing us to make dielectric tunability measurements. A counterelectrode beneath the thin film acts as a 
ground plane. 



III. PHYSICAL MODEL FOR THE SYSTEM 



A. Calculating the Electric Field Near the Probe Tip 



In order to arrive at quantitative results, we developed a physical model for the system, starting with the simplest 
case of a uniform bulk sample_J3ecause the probe tip length is much less than the wavelength A ~ 4 cm, we are able 
to perform a static calculationcJ of the microwave electric fields. Cylindrical symmetry further simplifies the problem 
to two dimensions. Because of the complicated geometry of the probe tip in contact with a multi-layered sample, we 
use a finite element model on a grid. We calculate the potential $ in the region represented by the grid, taking into 
account any changes in the permittivity e^, by solving Poisson's equation: 

V2$ + i-(v$) • (Ve^) = 0. (1) 

Using a rectangular grid, we represent the probe tip as a cone with a blunt end of radius tq (see Fig. H), and perform 
the calculation in a spreadsheet program. From Eq. (|l|), if $jj is the potential at the cell at column i and row j, 
is a weighted average of the values of the four adjacent cells: 

V A -\- ' 

r 

where r is the radius from the cylindrical coordinate axis, Ar and Az are the spacing between cells in the r and 
z directions, respectively, and we take Ar = Az. Equation (|^) is the simplest case; in practice, the equation is 
complicated by interfaces between the dielectric sample (e^ > 1) and the air (e^ — 1), for example. Near the probe 
tip (inside the dashed box in Fig. ||) where the electric field is the strongest, and hence, the most critical, the grid 
spacings Ar and Az are small and uniform. Inside this box we fix Az = Azi„ = 0.1/im, while the value of Ar = Ari„ 
is uniform, but variable (0.5 fim ^ Ar^ ^ l.O/xm), to allow for probes with different sharpnesses. We represent the 
probe sharpness with the aspect ratio parameter a = Azi„/Ari„. 

It is important for the boundaries of the grid to be sufficiently far away in order to minimize the effect of the chosen 
boundary conditions on the electric field near the probe tip. To accomplish this, outside the dashed box in Fig. ||, 
the values of Ar and Az continuously increase with distance away from the probe tip, allowing the outer radius of 
the grid to be at least 4 mm, and the height of the grid to be 2 mm. The resulting grid consists of 84 x 117 cells, 
which is small enough to be a manageable calculation with a modern personal computer. The top and outer boundary 
conditions were d^/dn = where n is the coordinate normal to the edge. At the bottom of the grid, which represents 
the bottom side of the 500 /xm thick sample, we used $ = for the boundary condition. To match this condition, we 
always place the sample on top of a metallic layer for scanning; this has the added benefit of shielding the microscope 
from the effects of whatever is beneath the sample, which could be difficult to model. 
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Two possible fitting parameters for the model are the aspect ratio a of the conical tip, and the radius rg of the 
blunt probe end. We obtain a satisfactory fit with our data (see below) by fixing tq = (0.6 fj,m.)/a. This leaves a as 
a single fitting parameter to represent all probes; typically we find that 1 < a < 2 for all of our probes. 

Shown in Fig. ||(a) is the calculated electric field near the probe tip for a probe with a = 1, and a sample with 
— 2.1. We notice that the fields are concentrated near the probe tip, as expected. The spatial resolution is related 
to the size of the probe tip, as shown by this concentration of the fields near the tip. Figure ^(b) and (c) show the 
calculated electric field as a function of radius away from the center of the probe, and depth in the sample, respectively, 
for two samples with = 2.1 and 305. For higher e^, the fields are more highly concentrated, and fall off more quickly 
away from the probe, indicating that the spatial resolution of the microscope is higher for high materials. From 
Fig. ||(b) and (c), the depth resolution appears to be slightly better than the lateral resolution. 



B. Calculating the Frequency Shift of the Microscope 

Using perturbation theoryj^^ we calculate the frequency shift of the microscope as a function of the fields near the 
probe tip. We define e^i and £^2 as the permittivities of two samples, the subscripts 1 and 2 indicating the unperturbed 



and perturbed system. If Ei and E2 are the calculated, 
the microscope upon going from sample 1 to sample 2 isl 



A/ _ 
f ^ AW 



ectric fields inside the two samples, the frequency shift of 
(e,2 - eri)Ei ■ E^dV. (3) 



Vs 



where W is the energy stored in the resonator, and the integral is over the volume Vs of the sample. We calculate 
an approximate W using the equation for the loaded Q of the resonator, Ql = ^oW/Pi, where luq is the resonant 
frequency, and Pi is the power loss in the resonator. In our case, we use a bare 500 /im thick LaAlOa (LAO) substrate 
for the unperturbed system (e^i = 24, -Bi),lj because its properties are well-characterized and it is a common substrate 
for oxide dielectric thin films. 

In the case of a bulk dielectric sample, we calculate the frequency shift A/ upon replacement of the LAO substrate 
by the sample of interest with permittivity er2- In the case of thin films on a LAO substrate, we calculate the frequency 
shift A/ associated with replacing a thin top layer of the LAO substrate with a thin film of permittivity 6^2 • 



IV. MEASURING THE PERMITTIVITY OF BULK SAMPLES 



To establish the validity of our microscope and model for performing permittivity measurements, we started with 
bulk dielectric materials. In order to obtain a comparison between the model presented aboV|&,and experimental 
results, we scanned a series of 500 /im thick bulk dielectrics with known microwave permittivities :EZI Teflon (e^ = 2.1), 
MgO {er = 10), LAO (e^ = 24), and SrTiOs (e^ — 300). In Fig. ^, the experimental frequency shift data points are 
shown for three different probe tips with different aspect ratios a. Model fits are also shown, where only the value of 
a has been varied. Note that the zero of frequency shift is for LAO (e^ = 24), the chosen unperturbed state of our 
resonator. 

To measure the permittivity of a dielectric sample, we must first calibrate the probe (i. e. determine the parameter 
a) using at least two samples with known and the same thickness. Once a has been determined from this calibration, 
we can measure of any sample with any thickness, using the finite element model and Eq. (^) to convert the measured 
A/ to Cr- 

We tested this method by measuring of a 500 laa thick KTaOa crystal, which is paraelectric with a cubic 
perovskitc structure at room temperature (T^ = 13 K)E5 WeJound the permittivity to be = 262 ±20, in agreement 
with microwave data in the literature (e^ = 240 at 9.4 GIIz)Ej and low-frequency data obtained with a parallel-plate 
capacitor (e^ = 260 at 100 Hz, and = 238.5 at 100 kHz). 



V. IMAGING PERMITTIVITY OF THIN FILMS 



For imaging thin films, we extended the finite element model to include a thin film on top of a dielectric substrate, 
as shown in Fig. |^. For the purpose of thin film imaging, we consider the unperturbed system to be the probe tip 
in contact with a bare 500 fim thick LAO substrate; the perturbed system, in addition, has a dielectric thin film on 
top of the 500 /im thick LAO substrate. Because the change in total thickness with the addition of the thin film is 
negligible compared to the 500 /im thick substrate, we can treat the thin film as the only perturbation to the system. 
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This perturbation due to the thin film is evident in Fig. ||, through the bending of the equipotential hues at the top 
and bottom surfaces of the thin film. Thus, the volume Vs in Eq. y includes only the volume of the thin film, and we 
calculate the frequency shift A/ associated with replacing a thin top layer of the LAO substrate with a thin film of 
permittivity £^2- 

The process for quantitative imaging of thin films is as follows. First, we determine the probe's a parameter using 
the method described above using bulk samples. Then, using the calculated electric fields from the finite element 
model for thin film samples, we use Eq. (|^) to calculate A/, in the same way that we did for bulk samples, except that 
we integrate only over the volume of the thin film. Using this model calculation, we obtain a functional relationship 
between A/ and 6^2 of the thin film, much like the results shown in Fig. ^ This relationship is used to convert A/ 
to er2, the permittivity of the thin film, after scanning the sample. 

To test the thin film model, we imaged thin-film samples of SrTiOs and Bao.eSro 4Ti03 (BST) on LAO substrates. 
The results are summarized in Table The thin-film permittivity was also measured at 10 kHz using Au interdigital 
electrodes deposited on the films. Both our microwave measurements and the int erdig ital electrode measurements 
are primarily sensitive to the in-plane component of e,. for these samples (see Sec. VI B| ). For the STO sample, the 
microwave permittivity is comparable to the low-frequency permittivity, showing that there is very little dispersion 
in this film. The BST filnas-jboth show significant dispersion, which is nonetheless within the range observed in the 
literature for similar films .EJO 



VI. OTHER MICROSCOPE ISSUES 



A. Spatial Resolution 

Using the finite element model, we can calculate the spatial resolution of the microscope. From Eq. (||) we see that 
the important quantity is the electric field dot product contained in the integral. Near the probe tip, the electric field 
is the strongest, and falls off nearly to zero at the outer boundary of the model grid. Thus, most of the contribution 
to the integral will come from the region in the sample near the probe tip. For our purposes, we define the lateral 
spatial resolution to be 2rres, where the integral over the volume V = nr'^^^d under the probe tip is equal to half of 
the integral over the volume of the whole sample. For thin films, the depth d is taken to be the thickness of the film, 
because we assume the substrate to be uniform and therefore not to contribute any features to the dielectric image. 
For bulk samples, we choose d = 2rres^ since we expect the vertical spatial resolution to be approximately equal to 
the lateral spatial resolution (see Fig. ||). 

Figure || shows the calculated spatial resolution for both a 500 /im thick bulk dielectric sample, and for a 400 nm 
thin film on a 500 /xm thick LAO substrate, for a typical probe with ro = 0.4 /im and an aspect ratio a = 1.5. We 
notice that the spatial resolution is the highest (2rres 1 to 1.3 /im) for thin films; this value for 2rres agrees with 
experimental results, and is illustrated graphically in Fig. ^ as a concentration of the electric field near the probe 
tip. At high permittivities, the spatial resolutions for bulk and thin film samples converge to about 1.5 /im, which is 
approximately twice the tip bluntness of 2ro = 0.8 /im. Model results show that the spatial resolution for any is 
approximately proportional to tq. 



B. Permittivity Tensor Directional Sensitivity 

Again using the model, we can determine the directional sensitivity of the microscope to the permittivity tensor by 
finding the relative magnitude of the radial (Ir) and vertical (/z) components of the integral in Eq. (^. This ratio is 
equal to 

Ir_ _ Ivs ^rl ■ Er2dV 

where Eri, Er2, Ezi, and Ez2 are the unperturbed and perturbed radial and vertical electric fields, respectively. This 
quantity Ir/Iz is shown in Fig. ^ as a function of sample permittivity for a 500 /im thick bulk dielectric and a 400 nm 
thin film on a LAO substrate. For low-permittivity bulk samples, the normal (z) component of the permittivity tensor 
dominates; for most thin films, the in-plane (radial) component dominates, as shown in Fig. ^, where the electric field 
in the thin film is mainly in the horizontal direction. 

We tested these model results by imaging a 5 x 5 x 10 mm 0.1% Cc-doped (Sra;Bai_2:)i-j/(Nb205)y (SBN) single 
crystal, with x = 0.61 and y = 0.4993. This anisotropic crystal has specified permittivities of en = 450 and €33 — 880.l3 
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We found that the microscope frequency shift was ~ 5% larger (more negative) when the probe was in contact with 
the face with €33 in-plane, relative to when the probe was in contact with the face with en in-plane; this agrees well 
with our model results, which predict a difference of 7%. 

C. Spherical Approximation for the Probe Tip 

As an alternative to our finite element model, an analytic solution for the electric field near the p rabe tip can be 
found by approximating the probe tip as a sphere of radius rsph, and assuming an infinite sample Jlj'ti We tested 
the usefulness of a spherical approximation for the probe tip by replacing the blunt cone with a sphere in our finite 
element model. To obtain the correct dependence of frequency shift on sample permittivity, it was necessary to include 
a second fitting parameter, the radius tq of a flat area of the sphere in contact with the sample. Figure ||(a) shows the 
best fit of A/ vs. er for one of our probes for 500 fim thick bulk samples using a spherical tip model with r^p/j = 14 /im 
and ro = 0.4 /im. 

The required radius of the sphere rsph is surprisingly large, so we calculated the spatial resolution of such a probe 
using the method described above, for a 400 nm thin film with Cr = 300 on a LAO substrate. Shown in Fig. ^(b) 
is the integral inside a cylinder of radius r^es as a fraction of the total integral in Eq. |[ for a spherical tip and a 
conical tip. For the sphere-shaped probe, the spatial resolution was 2rres = 3 /im. The conical probe model giving 
the best fit of A/ vs. has a = 1.3, giving a spatial resolution of 1.4 /im. In addition, because of the slow fall-off of 
the integral as a function of radius for a spherical tip [Fig. ||(b)], imaging with such a tip would result in the smallest 
image features being smeared out to a diameter > 10 /im. Since we do not observe this effect, and find our spatial 
resolution to he ^ 1.4 /im, we conclude that a cone is a much better approximation for the probe tip than a sphere, at 
least for the geometry of our system. We believe that this is due to the long-range nature of electromagnetism, where 
the potential due to a point charge diminishes,at the rate 1/r, causing the sides of the cone-shaped tip to have an 
important contribution to the frequency shift H Figure ||(a) illustrates this effect, by showing that even at a distance 
of 6 /im from the axis, the electric field is relatively strong at the cone surface. 

Also, a problem with an analytic solution using a sphere for the tip is the limitation to a single fitting parameter, 
the sphere radius rsph- For our probe tips, it was not possible to obtain a reasonable fit with the data on a A/ vs. 
er curve without an additional fitting parameter, such as the radius of a flat area on the sphere (a geometry which 
cannot be solved analytically). , ,■ , 

The assumption of an inflnite samplai3o can also pose problems. We find that this is an unrealistic assumption, 
because the properties of whatever is directly beneath the sample substrate can have a measurable effect on the 
frequency shift of the microscope, which will affect quantitative imaging. A better approach is to have a bulk 
conducting ground plane beneath the sample in the experiment and also in the model, so that everything beneath 
this ground plane can be safely ignored. 

VII. NONLINEAR DIELECTRIC IMAGING 

Electric field-dependent imaging can be accomplished by applying a voltage bias (VJ,) to the probe tip via a bias tee. 
in the resonator (see inset to Fig. |^). A metallic layer beneath the thin film acts as a grounded counterelectrode.E3 
In order to prevent the counterelectrade from dominating the microwave measurement (modeling of the system has 
shown this to be a potential problerro), we use a high-sheet- resistance counterelectrode, making it virtually invisible 
to the microwave fields and causing minimal frequency shift. This is demonstrated in Fig. 2(a) in Ref. |l^, where the 
frequency shift in the high-sheet-resistance limit is shown to saturate at a value dependent only on the substrate. As a 
result, the presence of the thin-film counterelectrode can be safely ignored in the finite element model described above. 
Because the counterelectrode is immediately beneath the dielectric thin film, the applied electric field is primarily in 
the vertical direction, unlike the microwave electric field, which is mainly in the horizontal direction for thin films 
with large permittivies (Fig. Also, by simulating the applied field using a finite element model similar to that 
presented above, we find that the applied electric field beneath the probe tip is approximately uniform and equal to 
Eb = Vb/tf, where tf is the thickness of the dielectric thin film. 

By modulating the bias voltage applied to the probe tip, we can extract nonlinear terms in the permittivity. 
Expanding the electric displacement D in powers of the electric field E, and keeping only the nonzero tcrmsJ13 we 
have 

Di{E) = eiiEi + -eii^EiEs + -en^^EiE^ + ... (5) 
where Ei is the rf electric field in the r direction, and E3 = Eb is the applied bias electric field in the z direction. 
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Adding a low-frequency oscillatory component {uti, = 1 kHz, with an amplitude of H = 1 V, in our case) taJthe 
bias voltage, the applied electric field is Ei, = Ef'^ + E^ cos uJbt. We find that the effective rf permittivity is thenla 



1 



(-rf — en 



reiia^^b'^ + eii33 
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+ (— H ^— 2-)£;{,cos(cjhi) 



1 

12 



(wizEl cos{2ujbt) 



(6) 



We note that the components of e^/ at lu{, and 2ujb are approximately proportional to ens and £1133, respectively. 
Expanding the resonant frequency of the microscope as a Taylor series about /o(er/ — en), we have 



/o[er/W] =/o(en) + 



dfo 



dCrf 



[(rf{t) - en] + 



(7) 
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Substituting Eq. ^ into Eq. 0, and keeping only the larger terms. 
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(8) 



Thus, the components of the frequency shift signal at ujb and 2ujb can be extracted to determine the nonlinear 
permittivity terms £113 and eii33. These nonlinear terms can be measured simultaneously with the linear permittivity 
(en) while scanning. 

As an alternative, the electric field could be applied in the horizontal direction using thin film electodes deposited 
on top of the dielectric thin film. The advantage in this case is that diagonal nonlinear permittivity tensor terms 
could be measured, such as em and enn; the disadvantage is that imaging is limited to the small gap between the 
electrodes. 



VIII. AN APPLICATION: MEASURING RECOVERY OF DIELECTRIC PROPERTIES IN A THIN 

FILM AFTER ANNEALING 

To demonstrate the usefulness of the microwave microscope, we scanned a sample consisting of a laser-ablated 
370 nm thick Bao.6Sro.4Ti03 (BST) thin film on a 70 nm Lao.95Sro.o5Co03 (LSCO) counterelectrode. The substrate 
is LaA103 (LAO). The films were deposited at 700 °C in 200 mTorr of O2. The sheet resistance of the LSCO 
counterelectrode is about 400 fi/D, sufficiently large to render it invisible at microwave frequencies. 

Figure ^(a) shows a schematic diagram of a 76.5 x 20 fim^ region of the thin film sample which we scanned. The 
gray-shaded areas indicate regions which were milled through the BST layer using a gallium focused ion beam (FIB). 
There is a 1 /im wide line and a corner of a 5 ^m wide "frame" surrounding an untouched 20 x 20 ^m^ region. Fig. 
^(b) shows a permittivity image of the region sketched in (a) . The narrow line shows up as a_wider band with low 
permittivity, a sign that the region near the line was damaged by the gallium ion beam tailso during milling. The 
wide frame appears as a double line in (b) due to edge effects at the edge of the milled area; also, the value of e^ 
shown in this wide frame region is invalid, since there is no dielectric thin_film there. Several low-permittivity regions 
appear randomly scattered over the scan area. We have shown previouslytj that these are particles ( "laser particles" ) 
on the surface which probably accrued during pulsed laser deposition. 

After acquiring the image shown in Fig. ^(b), we annealed the sample at 650 °C in air for 20 minutes. Afterward, 
we scanned the same region again; the resulting image is shown in (c). The overall thin film permittivity has increased 
between (b) and (c), and the narrow line has become less prominent, indicating that perhaps the permittivity of the 
damaged region was partially restored. 

Figure ^(d) and (e) show nonlinear dielectric images (eii3) taken with an applied dc bias of -3.5 V, giving an 
average vertical dc electric field of -95 kV/cm beneath the probe tip. Figure ||(d) was acquired before the sample was 
annealed, and shows that the tunability has been destroyed near the narrow FIB milled line, and within about 12 ^m 
of the FIB milled frame. Figure ^(e), which was acquired after annealing, shows that tunability has been restored in 
the damaged regions, and slightly improved in the undamaged regions. 
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Figure |l^ shows hysteresis loops taken after the sample was annealed, at the point marked by a "+" in Fig. 
^(c) and (e). The thin film permittivity is shown as a function of applied electric field in Fig. [lo|(a). As expected, 
the permittivity decreases when an electric field is applied. The small amount of hysteresis is probably due to the 
broadened transition temperature of the nominally paraelectric thin film. The first nonlinear dielectric term ens is 
shown in Fig. ^(b), while the second nonlinear term 61133 is shown in (c). Wc notice thaj_fii3 is close to zero for 
zero applied field, which would be expected for a material with no spontaneous polarization.Ej In addition, we expect 
the sign of 6113 to be sensitive to the direction of the polarization, which is evident in (b). The second nonlinear 
term, on the other hand, is nonzero at zero applied field. The curves in Fig. [l^ are centered at a nonzero field-c^ —20 
kV/cm, probably because the asymmetric capacitor electrodes induce unequal charges at the two electrodes.Elo The 
observed tunability of is srnail (~ 2%) probably because we are measuring an off-diagonal nonlinear component of 
the permittivity tensor (eii3).Ej 

We have evidence that the nonlinearity signal at Wf, (see Sec. VII) is related to the normal component of ferroelectric 
polarization, allowing us to image domain structures in bulk crystals. 



IX. SENSITIVITY OF THE MICROSCOPE 



We find the sensitivity of the microwave microscope by observing the noise in the dielectric permittivity and 
tunability data. For a 370 nm thin film on a 500 /im thick LAO substrate, with an averaging time of 40 ms, wc find 
that the permittivity sensitivity is Ae^ = 2 at = 500, and the nonlinearity sensitivity is Aeii3 = 10^'^ (kV/cm)~^. 



X. FUTURE IMPROVEMENTS 



The sensitivity of the microscope to sample permittivity could be improved by decreasing noise and drift in the 
electronics, using a high-stability microwave source, for example. A piezoelectric scanning system will reduce vibration 
noise and mechanical drift caused by stepper motors. The spatial resolution is currently limited by the probe tip 
radius, which is made worse by contacting the sample; it could be improved by operating the microscope out of 
contact with a distance feedback mechanism. The latter two improvements could reduce the line-to-line registration 
drift seen in Fig. ^. 
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TABLE I. Comparison of microwave permittivity Cr of thin-films on LaAlOa, measured with the microwave microscope, and 



low-frequency permittivity, measured using interdigital electrodes. 



Film material 


Thickness (nm) 


Microwave er 


Cr at 10 kHz 


SrTiOa 


440 


145 ± 26 at 7.2 GHz 


150 


Bao.6Sro.4Ti03 


300 


388 ± 14 at 8.1 GHz 


700 


Bao.6Sro.4Ti03 


400 


573 ±27 at 8.1 GHz 


1030 



FIG. 1. Schematic of the microwave microscope. The inset shows a close-up of the probe and a thin-film sample, along with 
the probe bias voltage, Vt- 



FIG. 2. The geometry of the finite element relaxation model (not to scale). The left border represents the cylindrical 
coordinate axis. Inside the dashed box, the grid cell spacings, Ar and Az are constant, while outside, they continuously 
increase toward the boundaries. 



FIG. 3. (a) The calculated electric field magnitude near the probe tip (in the region indicated by the dashed box in Fig. g) 
for a sample with permittivity tr = 2.1, and a probe with aspect ratio a = 1 and tip bluntness ro = 0.6 fj,m. The electric field 
magnitude as a function of radius (b) and depth in the sample (c) are shown for samples with er = 2.1 and 305. 



FIG. 4. Model results (lines) vs. data (symbols) with 500 fim thick bulk dielectrics for three different probe tips with different 
values of the aspect ratio a. Frequency shifts are relative to a LaAlOa (LAO) sample, with = 24. 
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FIG. 5. The electric field near a probe tip in contact with a thin film with = 300 on a substrate with ~ 24. The region 
shown is 6 X 6.8 /ira^ . The equipotential lines are shown; electric field lines are perpendicular to these lines. 

FIG. 6. Spatial resolution of the microscope (2rres) as a function of sample permittivity, for a probe tip with aspect ratio 
a = 1.5 and radius ro — 0.4 ^m (see Fig. |^. These model results are given for a 500 /im thick bulk sample, and a 400 nm thin 
film on top of a 500 fj,m bulk LAO substrate (tr = 24). 

FIG. 7. The directional sensitivity of the microscope to the sample's permittivity tensor. The ratio Ir/Iz gives the relative 
contribution of the radial (in-plane) to the vertical (normal) components of for a sample with isotropic er- Results are given 
for a 500 fim thick bulk sample, and a 400 nm thin film on top of a 500 /im bulk LAO substrate. 

FIG. 8. (a) Data vs. model results, using the approximation of a sphere for the probe tip, with 500 fim thick bulk dielectric 
samples. The spherical tip parameters are the sphere radius r^ph ~ 14 fim, and contact area radius ro = 0.4 fim. (b) The 
contribution to the electric field integral in Eq. ^ inside a radius rint as a fraction of the total integral. Results are given for 
the sphere-shaped probe tip, and a blunt cone-shaped tip (with a = 1.3), for a 400 nm thin film of permittivity er = 300 on 
a 500 ^m thick LAO substrate. The spatial resolution for the spherical tip is ~ 3 /im, while for the conical tip it is a more 
realistic ~ 1.4 pm. 

FIG. 9. Images of a 370 nm thick Bao.6Sro.4Ti03 (BST) thin film on a 70 nm Lao.gsSro.osCoOs (LSCO) counterelectrode on 
a LaAlOs (LAO) substrate. All images are of the same 76.5 x20 ^m^ region, (a) Schematic diagram of the focused ion beam 
(FIB) milled regions in the scan region, (b) Scan of the thin film permittivity before annealing, (c) Thin film er after annealing 
at 650 °C in air for 20 minutes, (d) Thin film tunability as shown in the nonlinear dielectric term ens, before annealing. The 
low-tunability regions were damaged during FIB milling, (e) Nonlinear dielectric (ens) image after annealing, showing that 
the tunability has been restored. 

FIG. 10. Hysteresis loops at the location marked with a "+" in Fig. ^ Permittivity (a), first nonlinear permittivity term 
£113 (b), and second nonlinear term £1133 (c) of the thin film as a function of the applied electric field E^'^. 
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